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Abstract 

Constraints from searches for squarks and gluinos at the LHC at -v/s=8 TeV are 
apphed to the parameter space of the NMSSM with universal squark/slepton and 
gaugino masses at the GUT scale, but allowing for non-universal soft Higgs mass 
parameters (the sNMSSM). We confine ourselves to regions of the parameter space 
compatible with a 125 GeV Higgs boson with diphoton signal rates at least as large as 
the Standard Model ones, and a dark matter candidate compatible with WMAP and 
XENONIOO constraints. Following the simulation of numerous points in the mo— Mx/2 
plane, we compare the constraints on the sNMSSM from 3-5 jets -|- missing Enp chan- 
nels as well as from multijet + missing Et channels with the corresponding cMSSM 
constraints. Due to the longer squark decay cascades, lower bounds on ^'^^ 
leviated by up to 50 GeV. For heavy squarks at large mg, the dominant constraints 
originate from multijet + missing Ex channels due to gluino decays via stop pairs. 



1 Introduction 



One of the most important tasks of the LHC - besides the quest for the Higgs boson - 
is the search for new elementary particles like those predicted in supersymmetric (SUSY) 
extensions of the Standard Model (SM). So far the search for such SUSY particles (sparti- 
cles) has not been successful; the absence of corresponding signal events can be interpreted 
as lower bounds on sparticle masses (see [THS] for recent ATLAS publications of results at 
y/s=S TeV, [7HTT] for recent CMS publications of results at y/s=S TeV, |T2] and the web 
pages [13] and [H] for summaries of searches for sparticles by the ATLAS and CMS collab- 
oration). Clearly these lower bounds on sparticle masses are not model independent, since 
they depend on the sparticle couplings and decay cascades, and hence on a large number 
of unknown parameters. 

The large number of unknown parameters of SUSY extensions of the SM is greatly 
reduced if one assumes universal soft SUSY breaking terms at the Grand Unification (GUT) 
scale, which is also theoretically appealing. Such models are denoted as "constrained" , and 
this is the case of the cMSSM (constrained Minimal Supersymmetric extension of the SM). 
Since the various sparticle masses and couplings are now strongly correlated, constrained 
models often serve as useful benchmark scenarios. The bounds on sparticle masses can then 
be represented as bounds in the mo — M1/2 plane, where mo denotes the universal squark 
and slepton masses and M1/2 the universal gaugino masses at the GUT scale. Frequently 
the cMSSM with tan /3 = 10 and Ao = is used to this end (tan /3 being the ratio of the 
two Higgs vevs (Hu) / {Ha), and Aq denoting the universal soft SUSY breaking trilinear 
couplings at the GUT scale). 

However, most of the considered parameter space of the cMSSM with tan /3 = 10 and 
y4o = is neither consistent with the observation of a SM like Higgs boson near 125 GeV 
[TSIITG] . nor with the dark matter relic density as determined by the WMAP experiment [T7] . 
In recent publications [T8l - l2^ . the LHC bounds on sparticle masses have been applied to 
the cMSSM (or variants thereof as the NUHM with non-universal soft Higgs masses at the 
GUT scale), but with parameters consistent with a SM like Higgs boson near 125 GeV, 
and/or dark matter consistent with WMAP bounds on the relic density and XENONIOO 
limits [25] on the dark matter direct detection cross section. Generally, one finds that the 
bounds on sparticle masses obtained within variants of the cMSSM or in the NUHM are 
similar to the cMSSM with tan /3 = 10 and Aq = 0. 

The MSSM is not the only possible supersymmetric extension of the SM. The simplest 
supersymmetric extension of the Standard Model with a scale invariant superpotential, i.e. 
where the soft SUSY breaking terms are the only dimensionful parameters, is the Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) [SB]. A supersymmetric Higgs mass 
term fi, as required in the MSSM, is generated dynamically by a vacuum expectation value 
(vev) of a gauge singlet (super-)field S, and is naturally of the order of the SUSY breaking 
scale. The attractive features of the MSSM are preserved, like a solution of the hierarchy 
problem, the unification of the running gauge coupling constants at a Grand Unification 
scale, and a dark matter candidate in the form of a stable lightest SUSY particle (LSP). 

The additional coupling A of the Higgs bosons in the NMSSM makes it much easier to 
accommodate a SM like Higgs boson near 125 GeV [27HSD] and alleviates the corresponding 
"little fine-tuning problem" of the cMSSM since hghter scalar top quarks (top squarks) are 



1 



possible [27H221EIHMl[Mll3SllMllliailSlllSlliai^ The fermionic component of the gauge 
singlet superfield S extends the neutralino sector of the MSSM, and this can lead to more 
complicated sparticle cascade decays [311E2H5S]- 

Like in the MSSM, one can consider constrained versions of the NMSSM with universal 
soft SUSY breaking terms at the GUT scale. A SM like Higgs boson near 125 GeV can 
easily be obtained within the semi-constrained NMSSM in which, similar to the NUHM, the 
soft SUSY breaking Higgs mass terms (and the trilinear couplings involving the singlet S) 
are allowed to deviate from the soft SUSY breaking terms involving squarks or sleptons. In 
several recent publications the parameter space of the semi-constrained NMSSM compatible 
with a SM like Higgs boson near 125 GeV (and possibly an enhanced diphoton signal rate) 
has been discussed [30l[35ll38ll45ll57] . The LHC constraints from negative squark and gluino 
searches had to be estimated or had been left aside, unless trivially satisfied due to very 
heavy squarks and/or gluinos. LHC constraints from the runs at y/s = 7 TeV, using the 
razor variables at CMS |58] , have been studied within a more restricted version of the semi- 
constrained NMSSM allowing only the soft SUSY breaking singlet mass term to deviate 
from mo in [59]. In this case the regions in the parameter space corresponding to large A, 
low rriQ and M1/2, which are interesting from the point of view of a 125 GeV Higgs with 
low fine-tuning, are not viable. 

Hence it becomes interesting and important to re-analyse the LHC constraints on sparti- 
cle masses within the semi-constrained NMSSM, which is the purpose of the present paper. 
We focus on the most constraining squark and gluino search channels analysed by the AT- 
LAS collaboration for the y/s = 8 TeV run: Searches for final states with jets and missing 
transverse momentum final states with large jet multiplicities [1], and with an isolated 
lepton [2]. 

We confine ourselves to phenomenologically acceptable regions of the sNMSSM param- 
eter space with a Higgs boson near 125 GeV, a diphoton signal rate near or above its SM 
value, a dark matter relic density in agreement with WMAP constraints, and a dark matter 
direct detection cross section compatible with XENONIOO constraints. The latter constrain 
the neutralino sector (the mass and the couplings of the lightest SUSY particle, the LSP) 
which has some impact on the sparticle decay cascades. 

In the NMSSM, sparticle decay cascades can differ from the MSSM for various reasons: 

• The higgsino-like neutralinos and charginos can be lighter than in most realistic sce- 
narios within the MSSM (with a Higgs boson near 125 GeV and a dark matter relic 
density compatible with WMAP constraints). The additional singlet-like neutralino 
(singlino) can mix with the MSSM like neutralinos, implying more complicate spar- 
ticle cascade decays [3111S2HSS]- 

• The top squarks are typically lighter than in realistic scenarios within the cMSSM 
(with a Higgs boson near 125 GeV), implying gluino/squark cascade decays via top 
squarks which lead to multijet events, reducing the missing transverse momentum 
and the average transverse momenta of jets. 

Hence it is not clear to which extent the bounds in the rriQ — M1/2 plane obtained by 
ATLAS for the cMSSM with tan/3 = 10 and Aq = are applicable to the semi-constrained 
NMSSM; the results of our present study allow to answer this question quantitatively. 
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In the next section, we briefly review the semi-constrained NMSSM and discuss our 
choice for points in the rriQ — M1/2 plane. In Section 3 we describe the tools used for the 
Monte Carlo studies and list the applied LHC constraints. In Section 4 we present the 
resulting bounds in the mo — M1/2 and Mgquark — ^giuino planes, discuss the origin of the 
differences in the bounds within the semi-constrained NMSSM with respect to the cMSSM, 
and summarise our conclusions. 



2 The semi-constrained NMSSM 

The NMSSM differs from the MSSM due to the presence of the gauge singlet superfield 
S. In the simplest realisation of the NMSSM, the fiHuHd Higgs mass term in the MSSM 
superpotential VFmssm is replaced by the coupling A of S* to Hu and Hd, and a self-coupling 
hiS^. Hence, in this version the superpotential Wnmssm is scale invariant, and given by: 

Wnmssm = ^SHu ■ Hd + —S + . . . , (1) 

where the dots denote the Yukawa couplings of Hu and Hd to the quarks and leptons as 
in the MSSM. Once the scalar component of S develops a vev s, the first term in VFnmssm 
generates an effective /x-term with 

/ieff = A s . (2) 

The soft SUSY breaking terms consist of mass terms for the Higgs bosons Hu-, Hd-, S, 
squarks = (mil, c^jl), "Uiij, (ij/j and sleptons ii = {i^n, (^il) and ej^j (where i = 1 ... 3 is a 
generation index): 

-Co = m^Ji/^p + m|Ji7rfp + m||^P + m||gi|2 + m||M4P + m||(i4r 

+m||/,|2 + m|Je4|2, (3) 

trilinear interactions involving the third generation squarks, sleptons and the Higgs fields 
(neglecting the Yukawa couplings of the first two generations): 

- £3 = (^htAt Q ■ Hu u^^R + /ifoAb Hd ■ Qd-iR + KA^ Hd ■ L 63% 

+ XA^ Hu-HdS+ ^^A^ S'^ + h.c. , (4) 

and mass terms for the gauginos B (bino), W"' (winos) and (gluinos): 
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M1BB + M2 W'Wa + M^ J2 G"Ga 



a=l a=l 



+ h.c. . (5) 



The neutral CP-even Higgs sector contains 3 states Hi, which are mixtures of the CP- 
even components of the superfields Hd and S. Their masses are described by a 3 x 3 
mass matrix M^^ij- The neutral CP-odd Higgs sector contains 2 physical states Ai, whose 
masses are described by a 2 x 2 mass matrix In the neutralino sector we have 5 
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states x^i which are mixtures of the bino 5, the neutral wino W'^, the neutral higgsinos 
from the superfields Hu and Hd-, and the singlino from the superfield S. Their masses are 
described by a 5 x 5 mass matrix A^^oj^. Expressions for the mass matrices - after Hu-, 
and S have developed vevs Vu, Vd and s, and including the dominant radiative corrections 
- can be found in [26] and will not be repeated here. 

As compared to two independent parameters in the Higgs sector of the MSSM at tree 
level (often chosen as tan /3 and M^), the Higgs sector of the NMSSM contains six parame- 
ters 

A , K , , A^, tan/3 = v^/vd , /icff ; (6) 

then the soft SUSY breaking mass terms for the Higgs bosons fn^Hui "^^d ^^"^ "^1 
determined implicitely by M^, tan/3 and /icg. 

In constrained versions of the NMSSM (as in the constrained MSSM) one assumes that 
the soft SUSY breaking terms involving gauginos, squarks and sleptons are universal at the 
GUT scale: 

Mi = M2 = M^ = Ml/2 , (7) 

2 2 2 2 2 2 /o\ 

"^g- = "^ui = = m^^ = = , (8) 

At = Ab = Ar^Ao. (9) 

In the semi-constrained NMSSM considered here, one allows the Higgs sector to play 
a special role: the Higgs soft mass terms rn'jj^, mfj and m| are allowed to differ from ttIq 
(and determined implicitely as noted above), and the trilinear couplings Ax, A^ can differ 
from Aq. Hence the complete parameter space is characterised by 

A , K , tan/3 , fi^s , Ax , A^ , Aq , M1/2 , mo , (10) 

where the latter five parameters are taken at the GUT scale. 

Subsequently we are interested in regions of the parameter space with large NMSSM- 
specific contributions to the SM-like Higgs mass, i.e. large values of A (and n) and low values 
of tan/3, which lead naturally to a SM-like Higgs boson H2 in the 125 GeV range p7H5U] . 
We impose constraints from LEP [60] on the lighter mostly singlet-like Higgs boson Hi, 
which still allow for a Hi mass below 114 GeV if its coupling to the Z boson is reduced. 
For H2 we require 124 GeV < Mh, < 127 GeV, a2liH2)/4li > 1 and aif (i72)/cri| ~ 1 
in order to comply with the observations at the LHC. 

We have implemented these constraints into a modified version of the public code 
NMSPEC [61] inside NMSSMTools [62l|63]. (In the Higgs sector we have used two-loop 
radiative corrections from [61], and for the top quark pole mass we have taken mtop = 
173.1 GeV.) The constraints from i?-physics are those of the version 3.2.0 of NMSSMTools, 
which are easily satisfied for the regime tan /3 < 3 relevant here. The dark matter relic den- 
sity and direct detection cross section of the LSP Xi (the lightest neutralino) are computed 
by MicrOmegas [^5H57] implemented in NMSSMTools. However, due to the low values of 
tan /3, the SUSY contribution to the anomalous magnetic moment of the muon Aa^ is not 
large enough to resolve the discrepancy between the SM and its measured value. 

Leaving aside Aa^, many regions in the space of the parameters in (ITOl) satisfy all the 
above conditions. Hence we proceed as follows: We start with a "lattice" in the mg — M1/2 
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plane, i.e. numerous fixed values for rriQ and Mi/2- For each fixed (mo, M1/2) we choose 
the remaining parameters such that not only the above phenomenological constraints are 
satisfied, but also in such a way that the lighter top squark mass and fics are relatively 
small. 

Light stop quarks and low fi^s minimise the fine-tuning [?rH29t[MH3M3M36|[39Hi2Ml^ 
HHHinilST]. On the other hand, present constraints from searches for these sparticles should 
be satisfied. First, since gluinos with masses below ~ 1 TeV are excluded, constraints from 
gluino mediated stop production OElEHlEn] turn out to be satisfied. Constraints from 
direct pair production of top squarks i [6| [70l471] require m^^ > 400 GeV for LSP masses of 
~ 80 GeV, as found below (note that slightly stronger bounds assume branching ratios and 
neutralino/chargino masses within simplified models which are not valid here; see [75] for 
proposals for search strategies for light stops within the general NMSSM). In addition, we 
require fi^s > 120 GeV so that the lighter chargino masses and chargino-neutralino mass 
splittings comply with present constraints. 

For each such point on a lattice in the mo — M1/2 plane, we perform Monte Carlo 
simulations (~ 10^ events), apply the cuts described in the next section, and compare the 
resulting signal event numbers to present constraints. This allows to identify viable regions 
(up to error bars) in the mg — M1/2 and Mgquark ~ ^gimno planes within the semi-constrained 
NMSSM. 

3 Monte Carlo simulation, search channels and verifi- 
cation 

For the calculation of the matrix elements we use MadGraph/MadEvent 5 [76], which 
includes Pythia 6.4 [77] for showering and hadronisation. Matching of the differential jet 
cross sections is performed according to the prescriptions in [TS]. The sparticle branching 
ratios are obtained with the help of the code NMSDECAY [Tg] (based on SDECAY [HD]), 
and are passed to Pythia. 

The output is given in StdHEP-format to the fast detector simulation Delphes |81j . 
Inside Delphes, the anti-k(t) jet reconstruction algorithm [82] is used, with the jet recon- 
struction performed by Fast Jet [53] . 

The sparticle (squark and gluino) production cross sections are obtained by Prospino 
at next-to-leading order (NLO) [8H - I86] . The resummation of soft gluon emission is taken 
into account in the form of a correction factor estimated from [87], and the theoretical 
uncertainties from scale and PDF choices are obtained from [87] [88]. 

To the output from Delphes we apply cuts on final states with jets and missing transverse 
momentum from searches for supersymmetry at y/s = 8 TeV with an integrated luminosity 
of 5.8 fb~^ by the ATLAS collaboration [1], which give at present the strongest constraints in 
the mo — Ml/2 plane in the cMSSM, as well as cuts on final states with large jet multiplicities 
from [1] and one isolated lepton from [2] which could, a priori, be relevant for the NMSSM. 

In Table 1 we summarise the cuts corresponding to the search channels which lead to 
the most stringent constraints in the mo — M1/2 plane in the semi-constrained NMSSM 
(depending on mo and M1/2), and the 95% confidence level (CL) upper limits (UL) on 



5 



Channel 


C-tight [4] 


D-tight [4] 


E-tight j4j 


9j55 [Ij 


8j80 [U 




4 


5 


6 


> 9 


> 8 




160 


160 


160 














4 


4 


PtUi) > 


130 


130 


130 


55 


80 


PT{j2---iN,J > 


60 


60 


60 


55 


80 


mefT(incl.)> 


1900 


1700 


1400 






£^r7^cfr(Nj)> 


0.25 


0.15 


0.15 






95% CL UL on Nse 


3.3 


6.0 


9.3 


5.4 


4.0 



Table 1: Cuts and 95% CL upper limits on the number Nse of signal events at y/s = 8 TeV 
and 5.8 fb~^ integrated luminosity for the search channels leading to the most stringent 
constraints in the mo — M1/2 plane in the cMSSM or the semi-constrained NMSSM. [E^^^^, 
Pt, Ht and meff in GeV.) 

the number A^^^ of signal events beyond the expected background in the corresponding 
channels, for an integrated luminosity of 5.8 fb~^. (More details on the event selections can 
be found in [HUj; bounds from the searches including one isolated lepton did not lead to 
stronger constraints.) Here mefr(Nj) is the scalar sum of the transverse momenta of E^^^^ 
together with the leading N jets, meff(incl.) the scalar sum of the transverse momenta of 
E^'^^^ together with all jets with px > 40 GeV, and Ht the scalar sum of the transverse 
momenta of all jets with pt > 40 GeV without E^^'^^. 

We first verified the validity of our simulations in the framework of the cMSSM with 
tan/3 = 10 and Aq = 0: For points in the mo — M1/2 plane along the 95% CL exclusion line 
in [1] we determined the number of signal events in all search channels, divided them by the 
corresponding upper limits given in [1] , and computed a ratio R from the most constraining 
search channel (giving the largest value for R). If our simulations would coincide exactly 
with those in [1], we would obtain R = 1. The resulting ratio R is shown as function of 
mo in Figure 1, where we also indicate the most constraining search channel by colours: 
black (full) for C-tight, blue (dashed) for D-tight and green (dotted) for E-tight. The most 
constraining search channels, depending on mo, coincide with the information given in [1]. 

We see that R deviates from 1 by up to ±30%, which we take as uncertainty in the 
number of signal events after cuts due to our simulation (it is considerably larger than our 
statistical error). This error is added linearly to the error on production cross sections from 
scale and PDF choices (which are slightly squark and gluino mass dependent). For a given 
value of mo we determine three different values of M1/2: one such that the number of signal 
events in the sNMSSM coincides with the 95% CL UL in Table 1 (for the most constraining 
channel), and two more such that the number of signal events in the sNMSSM coincides 
with the 95% CL UL ± the relative error obtained as before. This leads to an exclusion 
curve in the mo — M1/2 plane, as well as to curves which represent our errors on the number 
of signal events. 
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Figure 1: The ratio R of the number of signal events A^^^; from our simulation within the 
cMSSM along the 95% CL exclusion line in [1], divided by the upper limits given in Table 1 
from [Ij. The colours indicate the most constraining search channel: black (full) for C-tight, 
blue (dashed) for D-tight and green (dotted) for E-tight. 



4 Results and discussion 

After simulation of a variety of points in the mo — M1/2 plane in the semi-constrained 
NMSSM and application of the cuts, we require that accepted points give a number of 
signal events below the 95% CL upper limit in each search channel, ± the uncertainty as 
determined above. This procedure generates the bounds in the mo — M1/2 plane shown in 
Fig. 2, together with the uncertainties indicated by dashed lines. (The colours in Fig. 2 
indicate the most constraining search channels: red for the jets + E^^^^ channels C-tight, 
D-tight or E-tight from [4], and blue for the multijet -|- E^'^^^ channels 9j55 or 8j80 from [1].) 

Note that, for a given value of mo, the upper and lower error lines can originate from 
different search channels. For comparison, we also show the bounds obtained by ATLAS 
for the cMSSM with tan /3 = 10 and Ao = as a black line. 

We see that for lower values of mg, the cMSSM bounds are alleviated due to NMSSM 
specific sparticle decay cascades. For larger values of mo the bounds within the semi- 
constrained NMSSM seem stronger. This is due to the fact that in this regime bounds from 
the multijet channels [1] become stronger than the bounds from the otherwise dominant 
channel E-tight; notice that the bounds from the multijet channels [1] are not included 
in the ATLAS bounds for the cMSSM with tan/3 = 10 and Aq = Q from g]. We have 
verified that the bounds from the multijet channels would also dominate in the cMSSM 
for mo > 1400 GeV, leading to cMSSM bounds somewhat stronger than those given in [1]. 

The corresponding bounds in the Mgquark — ^giuino plane are shown in Fig. [3l Here 
the region mo < 1500 GeV (M1/2 ^ 450 GeV), where the bounds from the jets + 
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Figure 2: Bounds in the mo — M1/2 plane in the semi-constrained NMSSM (dashed lines 
indicate our error bars), and ATLAS bounds for the cMSSM with tan/3 = 10 and Aq = 
from |1] in black. The colours indicate the most constraining search channel for the NMSSM: 
Red for the jets + E^^^^ channels C-tight, D-tight or E-tight from jl], and blue for the 
multijet + EJp^^^ channels 9j55 or 8j80 from [1]. 

E^^^^ channels are somewhat weaker in the sNMSSM than in the cMSSM, corresponds 
to M„k < 1700 GeV (M^m^o > 1200 GeV). 

The different cMSSM and sNMSSM bounds are due to having distinct sparticle decay 
cascades. In the following, we discuss these cascades (which depend on mo) in more detail. 

For mo < 1000 GeV, the sparticle production cross section is dominated by up/down 
squark pair production. In the cMSSM with tan/3 = 10 and Aq = 0, the /i-parameter 
and hence the higgsino masses are relatively large (~ 800 GeV), the binos and winos are 
approximately eigenstates with masses ~ 0.4 x M1/2, ~ 0.8 x M1/2, respectively, and the 
bino is the LSP (violating generally WMAP bounds on the relic density). The dominant 
decays of the right-handed squarks and left-handed squarks ql are 

QR^q + Xi; QL^q + xt ^Q + W^ + Xi, (11) 

where Xi is essentially bino-like and Xi" essentially wino-like. Hence the dominant decay 
cascades are relatively short. 

In the semi-constrained NMSSM with a SM-like Higgs mass of ~ 125 GeV and a dark 
matter relic density consistent with WMAP constraints, the effective /i-parameter (and 
hence the higgsino masses), as well as the singlino mass parameter 2ks, are relatively small, 
in the 115 — 250 GeV range. Apart from alleviating the "little fine-tuning problem", such 
higgsino and singlino mass parameters generate large mixing angles in the neutralino sector. 
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Figure 3: Bounds in the Mgquark — ^giuino plane in the semi-constrained NMSSM (dashed 
hues indicate our error bars), and ATLAS bounds for the cMSSM with tan /3 = 10 and 
Aq = from in black. Colour/line code as in Fig. [2l 



The LSP, the lightest neutralino, is a mixture of higgsinos and singlino. The mostly bino- 
like neutralino is i-e. not the LSP, and the lighter chargino xf is essentially higgsino-like. 
Hence the dominant decays of the right-handed and left-handed squarks are 

qR ^ q + Xl^ q + + ^ q + + W^* + X'l ; (12) 
qL^q + xf -^q + Z/H + x^ -^q + Z/H + W^* + xl (13) 

The squark decay cascades lead to considerably more final states in the NMSSM, implying 
less missing transverse momentum and less pt per jet compared to the cMSSM. This ex- 
plains the lower number of signal events, and the somewhat lower bound in the mo — M1/2 
plane. An example for a benchmark point with these properties (with mo = 600 GeV and 
Ml/2 = 650 GeV) is given in Table H 

For mo close to 250 GeV (and M1/2 ^ 800 GeV), another phenomenon appears in 
the semi-constrained NMSSM with a SM like Higgs mass of ~ 125 GeV: The Higgs mass 
requires a non- universal soft Higgs mass term rriH^, which is considerably larger than mg 
at the GUT scale. This has some impact on the running squark and notably on the slepton 
masses from the GUT to the weak scale, leading to light sleptons. In fact, the LEP2 
bound on light slepton masses of ~ 100 GeV leads to a lower bound mo > 250 GeV on 
the parameter space. Moreover, the mostly bino-like neutralino decays dominantly into 
sleptons in this region. We have checked that constraints from searches in the channels 
including isolated leptons in [2] are satisfied in this region, and the dominant constraints 
still originate from the channel C-tight. 

For mo > 1000 GeV, the sparticle production cross section becomes dominated by 
squark-gluino production and, for mo > 1400 GeV, by gluino pair production. In the 
cMSSM, gluinos undergo three-body decays involving virtual squarks, with some preference 
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A (Msusy) 


0.64 


Mh, 


123 


(A^'susy) 


0.36 


Mh, 


125 


tan/3 (Msusy) 


3.02 


Mh, 


377 


/ioff (Msusy) 


120 


RT (ggF) 


1.73 


Ml/2 


650 


Rt (ggF) 


1.05 


niQ 


600 


M^o (LSP) 


76.5 


Ao 


-1262 


Hd comp. or Xi 


0.49 


A 


-4zo 


Hu comp. of Xi 






-176 


S comp. of x!i 


0.43 


{ "^^squarks 


1425 




159 


-^^gluino 


1490 




196 


Mi, 


604 


M^o (bino) 


283 




0.0941 


M^± (higgsino) 


112 


p 


8.4 X 10"^° 


M^o,M+ (winos) 


540 



1 IQ/^ m ro 


Ijrti^ /o J 


Ml Xs + w 


31 


ul X2 + " 


«o 


"L — ^ X5 + " 


oz 


di-^ X2 +u 


00 


ur Xa + u 




J V 1 J 

"i? X4 + « 


94 


(1 H- TT7-ZIZ 

X4 xf + 


69 


, 1 rr 


18 




45 


, 1 ^7 

xt^X2 + z 


17 


, 1 rr 


15 




100 


X? ^ xf + ^ 


23 


x^ ^ xf + ^1 


20 


xr ^ x^ + 


21 


xr -> x^ + 


17 



Table 2: Input parameters, spectrum and some branching fractions of a benchmark point 
with rriQ = 600 GeV, M1/2 = 650 GeV. All dimensionful quantities are given in GeV. a^j 
denotes the spin independent LSP-proton cross section, for which the present XENONIOO 
bound is < 3 X 10~^ for a LSP mass of ^ 77 GeV. Rj^ , -Rf ^ denote the signal cross sections 
of H2 relative to the SM, which are given in the gluon fusion production mode (ggF). The 
reduced signal cross sections of Hi in these channels are about 0.25. Hu,d and 5* denote the 
higgsino and the singlino components of Xi? respectively. 

for the somewhat lighter squarks of the third generation. In the semi-constrained NMSSM, 
the non- universal soft Higgs mass term rriH^, as well as the larger value for the top Yukawa 
coupling ht (due to the lower value of tan/3), lead to lighter stop masses. Apart from 
alleviating the "little fine-tuning problem", such stop masses imply dominant gluino two- 
body decays into the top quark + the lighter top squark (practically 100%). The latter 
decays dominantly into a bottom quark + the lighter chargino. (As discussed above, we 
verify that present bounds from top squark searches are satisfied.) Hence the final states 
involve a large number of jets, but a somewhat reduced missing transverse momentum. 
Whereas many of these final states would pass the cuts in the searches for jets and missing 
transverse momentum in [1], the searches for large jet multiplicities in [1] become now more 
relevant and, in fact, the constraints from the multijet channel are now dominant. For this 
reason, the constraints for tuq > 1400 GeV in the semi-constrained NMSSM are stronger 
than those from the channels D-tight and E-tight on the cMSSM in [1]. This would not be 
the case if the constraints from the multijet channel would be applied to the cMSSM. A 
second benchmark point with these properties (with mo = 1900 GeV and M1/2 = 450 GeV) 
is given in Table [31 
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A (Msusy) 


0.61 




95 


(Msusy) 


0.37 




126 


tan/3 (Msusy) 


2.31 




331 


^ i\/r \ 

f^es (Msusy) 


lib 


Ri (ggF) 


1.78 


1 /2 


^ou 


^2 ys&^ ) 


1 An 

1 .^u 


''''U 


1900 


M fLSP) 


79 


An 


-875 


; rrnnTi of 'v'^ 

11^ L,Ullip. Ui Xi 


54 




-296 


Hy^ conip. of ^^j* 


0.68 




-385 


S comp. of Xi 


0.32 




2070 




162 


-'^gluino 


1140 




192 




555 


M^o (bino) 


212 




0.105 


M^± (higgsino) 


109 


P 
^S/ 


7.2 X 10-9 


M^o,M+ (winos) 


398 



Decays 


T3T) {0/\ 

DtX[ /o ) 


g ^h+t 


100 


h^xt + b 


51 


ti^Xi + t 


22 


h ^X2 + t 


18 


^ X? + z* 


79 


xf + w* 


20 


xf -> + w* 


100 



Table 3: Input parameters, spectrum and some branching fractions of a benchmark point 
with niQ = 1900 GeV, M1/2 = 450 GeV. The reduced signal cross sections of Hi in the 
channels bb, 77 and ZZ (in the gluon fusion production mode) are ~ 0.025. 

We conclude that the present bounds on mo and M1/2 in the semi-constrained NMSSM 
with a SM like Higgs mass of ~ 125 GeV are somewhat alleviated with respect to the 
cMSSM for mo < 1000 GeV (mostly for mo ~ 500 GeV) due to the longer and more 
complicated sparticle decay cascades. This phenomenon could have been anticipated. Here 
we have studied it concretely with the result that the lower bound on M1/2, at fixed mo, 
decreases by at most 50 GeV for mo ~ 500 GeV; for mo > 1100 GeV (Msq^ark > 1500 GeV), 
the lower bound arises from multijet searches, whereas those of the cMSSM were derived 
from jet and missing transverse momentum search channels. The central line in Fig. 2 
serves to compare the bounds to the cMSSM whereas, to be conservative, the lower dashed 
line including our uncertainties should be used for bounds in the mo — M1/2 plane in the 
semi-constrained NMSSM. Moreover it is likely that, within the general NMSSM (as in 
the phenomenological MSSM in [23]), lower bounds on sparticle masses are much weaker. 
We expect that a large variety of corresponding scenarios are possible within the general 
NMSSM, which will require more dedicated studies. 
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